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Abstract. This paper demonstrates that random, independently chosen equi-dimension- 
al subspaces with a unitarily invariant distribution in a real Hilbert space provide nearly 
tight, nearly equiangular fusion frames. The angle between a pair of subspaces is measured 
in terms of the Hilbert-Schmidt inner product of the corresponding orthogonal projections. 
If the subspaces are selected at random, then a measure concentration argument shows that 
these inner products concentrate near an average value. Overwhelming success probability 
for near tightness and equiangularity is guaranteed if the dimension of the subspaces is 
sufficiently small compared to that of the Hilbert space and if the dimension of the Hilbert 
space is small compared to the sum of all subspace dimensions. 



1. Introduction 

A collection of closed subspaces in a Hilbert space is a fusion frame if a weighted sum 
of the corresponding orthogonal projections provides an approximate identity. Research 
on fusion frames has enjoyed a rapid growth in the frame theory literature of the last 
decade, see |12| and references therein. There are many applications of this field, driven 
by demands from distributed sensing [9], parallel processing [3], communication theory 
[1], quantum computing [6\ and even neuroscience |34j . Fusion frames are important for 
these applications because they model linear, distributed signal representation strategies. 
A primary goal in many of these settings is to make the representation robust to partial 
data loss when a signal to be stored or transmitted is projected onto the subspaces given by 
a fusion frame. If the components in the subspaces are lost then this constitutes an erasure. 
The performance of a fusion frame is typically measured by its ability to compensate for 
such lost components, either in a deterministic, adversarial erasure regime or in an averaged 
sense with additional statistical assumptions on the signal and on the erasures. Earlier work 
on fusion frame design and erasures has shown that typical notions of optimality have a 
geometric characterization: If equi-isoclinic fusion frames exist, then they are optimal for 
the worst-case error for recovery based on the canonical dual when up to two subspace 
components are removed, and if equi-distant fusion frames exist, then they are optimal for 
the mean-squared error after applying a Wiener filter to suppress the effect of erasures, see 
also a simpler setting in which they provide optimality [5j. The existence and construction 
of such subspace collections depends on the dimensions and has remained a challenge despite 
many efforts [32l [29l [28l [23l [25l [26l [271 El [II] . The motivation for this paper is to show that 
a random choice of subspaces is with high probability nearly tight and nearly equiangular. 
In a certain regime, this was already accomplished in an earlier paper by the construction of 
nearly tight frames with nearly isoclinic subspaces from matrices with the restricted isometry 
[7], but this only resulted in nearly orthogonal subspaces. It somehow seems natural that 
a random selection favors mutual orthogonality which is in a way the simplest case of 
equiangularity. A generalized form of the Welch bound shows that an equiangular, equi- 
dimensional tight fusion frame minimizes the maximum value for the Hilbert-Schmidt inner 
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product of any pair of orthogonal projections P and P' corresponding to two subspaces, 
resulting in the value of the inner product [311 Q] 

As N, K,s ^ oo, A r tr[PP']/s 2 — > 1. We establish these asymptotics for a random choice of 
subspaces. We cover this case of non-orthogonal equiangular subspaces by further develop- 
ing the partial orthonormalization strategy applied to frames that can be partitioned into 
nearly tight Riesz sequences from [7j. Such frames are implicitly constructed in the com- 
pressed sensing literature as the column vectors of sensing matrices that have the restricted 
isometry property. Many of the standard construction methods rely on randomization, 
meaning they pick a random sensing matrix which is shown to have the restricted isome- 
try property with high probability. In the present paper, we use a similar strategy with a 
specific choice of random frames and investigate the properties of the corresponding fusion 
frames, given by unitarily invariant, independently selected subspaces. We conclude that 
these are with overwhelming probability nearly tight and nearly equidistant, and hence 
near-optimal for many applications. 

The main results are summarized as follows, see the next section for notation: 
If {Vj}j" =1 are s-dimensional subspaces that are chosen independently at random in a 
iV-dimensional Hilbert space, with a unitarily invariant distribution, and e, 5 > 0, 1 + e = 
(1 + <5) 6 , then with probability at least 

1 _ 2(1 + - \N e -MS*/4+M5 3 /3 _ 2K( , + f y e -N5y4+N53/3 

5 S 

they form an e-nearly Ks /M-tight fusion frame. 

Under the same assumptions, letting each Pj denote the orthogonal projection onto Vj, 
and with 1 + e = (1 + 5) 3 , the subspaces form a nearly equiangular fusion frame, meaning 
for all j ^ I, 



1 (l + e)N N . , , /(l + e)AT Ne 2 



1 + eV s s 2 V s 4s 

where the failure probability is bounded by the sum of 

K(K + 1) se (l+<5W2-s(s-l)(5 2 /2-5 3 /3)/2 

and 

K(K + 1)((M + K(l + l)^) e -^ 2 /4+^ /3 + (1 + 4 )We _ M52/4+M5 3 /3) 



Examining the terms in these error bounds shows that as N, s, K — > oo and s/N — > c > 
then the failure probability decays exponentially in N if 

^ + _i M1+ 4 )<4V4 _ aV 3 

and 

^ln(l + ^)<5 2 /4-<5 3 /3. 

This paper is organized as follows: In the next section, we fix notation and recall ele- 
mentary results. Section [3] demonstrates that random subspaces lead to nearly tight fusion 
frames. The final section is dedicated to showing that the subspaces are nearly equiangular. 



random fusion frames are nearly equiangular and tight 3 

2. Frames, Fusion Frames and Riesz sequences 

We briefly review frames, fusion frames and Riesz sequences. 

Definition 2.1. A family of vectors {<Pj}i£j in a real or complex Hilbert space % is a frame 
for H if there are constants 0<A<B<ooso that for all x £ % we have 

A\\x\\ 2 <^2\{x,ipj)\ 2 < B\\x\\ 2 . 

If A = B then we say that the frame is A-tight, and if A = B = 1, it is a Parseval frame. 
If there is c > and \\<Pj\\ = c for all j 6 J then it is called an equal norm frame, and if 
c = 1 it is a unit norm frame. The analysis operator of the frame T : % — > £%(J) is given 
by (Tx)j = {x,(fj). The synthesis operator T* is the (Hilbert) adjoint of T and the frame 
operator is the positive self-adjoint invertible operator S = T*T. 

We recall that if {tpj}j^j is a frame with frame operator S then {<S , ~ 1 ' 2 ¥>i}ie.7 is a Parseval 
frame for rl. 

While frames assign scalar coefficients to a vector, fusion frames map it to its components 
in subspaces [17] . 

Definition 2.2. Given a real or complex Hilbert space % and a family of closed subspaces 
{Wi}i£K with associated positive weights < Vi, i £ K, then {Wi,Vi}i^x is a fusion frame 
for H if there exist constants < A < B < oo such that 

O ^ ^ Oil 1 1 1 1 

^4|MI < / y Vj \\Pi%\\ < B\\x\\ foranyxGH, 

i&K 

where each Pi is the orthogonal projection onto Wj. A fusion frame is called tight if A and 
B can be chosen to be equal, and Parseval if A = B = 1. For e > 0, the fusion frame is 
e-nearly tight if there is a constant C so that A = j^C, B = (1 + e)C. The fusion frame 
is equi-dimensional if all its subspaces Wi have the same dimension. If {Wi}i & x are closed 
subspaces ofH, we define the space 

0VFi = {ip = (ipi) ieK \ A G Wi,(ip,ip) < oo}, 
i€K 

with the inner product given by 

ieK 

The analysis operator of the fusion frame is the operator 

T-.n^^Wi, 

given by 

Tx = (viPix) ieK . 
The synthesis operator of the fusion frame is T* and is given by 

T*(ipi) i( z K = Viipi. 
ieK 

In analogy with frames, the fusion frame operator is the positive, self-adjoint and invertible 
operator S = T*T. 

We also need the notion of e- Riesz sequences. We choose the convention from [7J for tight 
Riesz sequences which is convenient for orthonormalization. 
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Definition 2.3. A family of vectors {}fi}f = i in a Hilbert space % is a Riesz basic sequence 
with lower (resp. upper) Riesz bounds < A < B < oo if for all scalars {a{\f =l we have 

N N N 

i=l i=l i=l 

This family of vectors is e-Riesz basic if for all scalars {ai}^ =1 we have 

^ N N N 

-^~ e l a *l 2 - II X^^ll 2 - ( x + e ) l ai l 2 - 

i=l i=l i=l 

3. Random fusion frames are nearly tight 

We first recall that random Gaussian vectors form a frame that can be partitioned into 
nearly orthonormal systems. 

Lemma 3.1. Let X be a random N x M matrix whose entries are independent, standard- 
normal distributed random variables, and let Xj denote the vector containing the entries of 
the jth row. Let u G M. M , \\u\\ = 1, and Z = Yl!j=i \( u ^j)\ 2 then 



J :r- 

-A r <5 2 /4+A f 5 3 /6 



P(Z > 1 + 5) < e 
and 

P(Z < — ) < e -^ 2 /4+N6V3 

1 + 6 

Proof. The distribution of the row vectors is unitarily invariant, so without loss of generality, 
we can set u = e±, the first vector of the canonical basis for M. M . The sum Z = Ylf=i -^f l 
is up to the normalization factor 1/N chi-squared distributed with N degrees of freedom. 
The usual combination of the Laplace transform and the Chernoff bound gives measure 
concentration. We have 

nz > i + s) < (i - t )-m e -tN(i+sy2 

for any < t < 1, so after choosing t = 5/(1 + 5) and truncating the Taylor expansion of 
the exponential 

P(Z > 1 + 5) < e ^Wl+5)/2-iV<5/2 < e -N5y4+N63/6 _ 



Similarly, for t > 0, 



i( Z < _J_) < e -Nln{l+t)/2+tN(l+S)- 1 /2 



so setting t = 5 and comparing the terms of the Taylor series in the exponent gives the 
desired bound. □ 

Choosing u among the M canonical basis vectors in M M together with a union bound 
shows that the norms of all the columns are nearly a constant. 

Corollary 3.2. With the random N x M matrix X as above, 

1 1 N 

P(^— <ly X ^<l + 5 for all I) > 1 - 2Me- N52 / A+N53 / 3 . 
1 + 5 ~ N ^ JJ ~ ' 
i=i 

We derive a stronger consequence with an argument similar to the exposition in Baraniuk 
et al. [2], as presented in [5]: With high probability, sufficiently small subsets of the column 
vectors of X have nearly tight upper and lower Riesz bounds. 
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Lemma 3.3. Let W = spanje^ , Ej 2 , . . . Ej s } in M. M , let {Xj}^ =1 be a random family of 
vectors in ]R M as above and let < 5 < 2, then the set 

X<={X: Z(y) < |M| 2 for all y G W} 

defined by the random variables Z(y) = J2f=i \ {Vi^j)\ 2 f or U £W has probability 

P(A ' < )< ( l + l re -^ 2 /4+^3 / 3_ 


Moreover, if 5 < 1, then 

X> = {X : Z{y) > (1 + 5f\\yf for all y G W} 
has the same upper bound for its probability as X< . 

Proof. Using the triangle inequality and Lipschitz continuity of the norm we can bootstrap 

from a net S with min ||y — < - for all y G W, \\y\\ = 1. By a volume inequality for 

toes' 2 

sphere packings, we know there is such an S with cardinality 

ws(i+f)'. 

Applying a union bound for the probability of the complement of 

X s = {X : Z(w) > — — Alwll 2 for sdlw G S} 

we get 

¥(x s ) > l - (l + ±y e - NS2 /*+ N53 / 3 

o 

Now let a be the smallest number such that Z 1 l 2 (y) > WvW holds for all y G W. We 

show 1 + a < (1 + <5) 3 . To see this, let y G W, \\y\\ = 1 and pick w G S, \\x — w\\ < |. Then, 
using the triangle inequality yields 

z v 2w > |z v V) _ z v 2to _ m)l >^__^_. 

Since the right hand side of the inequality chain is independent of y, according to the 
definition of a we obtain 

1 1 5/2 

(l + a)V2 " (1 + 5)1/2 " (l + a) 1 / 2 ' 

Solving for (1 + a) -1 and further estimating gives 

1 1 
T+a~ ~ (l + <5) 3 ' 
For the second inequality, choose 

X' s = {X : Z{w) < (1 + <5)||u;|| 2 for all it; G S} 
and establish the same bound for its probability as for Xs- Let a be smallest such that 

Z 1/2 (y) < {1 + a) l / 2 \\y\\ ,y G W. 
We again use the triangle inequality to obtain that if X G X' s then for any y G W with 

IMI = i 
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Again by definition 



(l + a) i/2< ( i + 5)1/2 + (1 + a) i/2* 



so 

l + <5 

l + a< 



1 - t5/2) 2 

and if 5 < 1 then (1 — 5/2)~ l < 1 + S, which implies 

1 + a < (1 + 5) 3 . 

□ 

A union bound gives a lower bound for the probability that all subsets in a partition of 
the column vectors of X have good Riesz bounds. 

Corollary 3.4. Let X be a random N xM matrix whose entries are independent, standard- 
normal distributed random variables, < 5 < 1, and let {1,2, ...,M} be partitioned into 
sets {Jk\k=i of maximal size max^ \Jk\ < s, then with probability 

l_2K(l + l) s e -^ 2 /4+^/3 > 


for each k the set of resettled column vectors {(Xjj/\nSl)jL 1 }i e j k forms an e-Riesz sequence 
with e = (1 + 5 f - 1. 

Proof. If the matrix is multiplied by the normalization factor, then a fixed set of columns 



{(Xjj/\/N)jL 1 }i£j k indexed by Jk in the partition with | J^\ < s is by the preceding lemma 
e-Riesz with 1 + e = (l+c5) 3 with probability of at least 1-2(1 + |)s e -JV 1 5 2 /4+Jva 3 /3_ Applying 
the union bound for all K sets in the partition gives the claimed estimate. □ 

The same proof allows us to establish frame bounds if we think of X as the analysis 
operator. If we choose a trivial partition with s = M < N then the above lemma states that 
the family of row vectors {Xj}^ =1 with standard normal entries forms with high probability 
a nearly tight frame for M M . We swap iV and M for later notational convenience. 

Corollary 3.5. Let X be a random MxN matrix whose entries are independent, standard- 
normal distributed random variables, then with probability at least 

1 _ 2(1 + l)^ e -M5 2 /4+M5 3 /3 _ 

5 

the row vectors {Xj}^L l of X form a frame for E> N with lower and upper frame bounds 

M M 
j=l 

Theorem 3.6. Given a family of random, independent subspaces {Vk}f =1 of M. N , whose 
distribution is invariant under unitaries and whose dimensions Sk = dim are bounded by 
Sk < s for all k, and let M = Y2k=i s k then with probability at least 

1 _ 2(1 + ^\N e -MS 2 /4+M6 3 /3 _ 2K <1 + ^y e -N8 2 /A+N8 3 /3 

5 5 

the fusion frame {Vj, l}JLi has upper and lower frame bounds M/(N(1 + 5) 6 ) and M(l + 
5f/N. 
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Proof. The preceding lemma implies that the random Gaussian M x N matrix is 
with high probability the analysis operator of a frame that is nearly M/iV-tight, because 

1 M 

M/(N(1 + 5f)\\x\\ 2 < - |(x,X,)| 2 < M(l + 5f/N 

i=i 

for all x£R N except for a set of probability 2(1 + |)JV e -M5 2 /4+M<5 3 /3. 

Moreover, let {1, 2, ... , M} be partitioned into subsets { Jk}k=i °f s i ze rnaxj |Jj| < s, then 
a union bound shows that with probability bounded below by 

l_2K(i + i)* e -^ 2 /4+^3/3 



each set of rescaled row vectors {Xj/V^V}jej fe has upper and lower Riesz bounds (1 + 5) 3 
and (1 + 5) -3 , respectively. 

Orthonormalizing the Riesz sequences then changes the frame bounds by at most a 
factor of (1 + <5) ±3 . The resulting frame is by construction partitioned into orthonormal 
systems, and thus equivalent to a fusion frame with the same frame bounds. Combining the 
probabilities for the failure of the frame bounds and of the Riesz bounds gives the stated 
bound. □ 



4. Random fusion frames are nearly equiangular 

Using a variation of the strategy by Dasgupta and Gupta [23], we obtain measure con- 
centration for an average of projected norms, which implies that with a proper scaling of 
dimensions, random subspaces become nearly equiangular. 

Lemma 4.1. Given {Xjj : 1 < j < s, 1 < I < N}, independent identically standard-normal 
distributed random variables, N > s, and < /3 < 1 then 

s s N 

i=i i=i i=i 
Similarly, if (3 > 1, then 

s s N 

j=l 1=1 1=1 

Proof. The first probability under consideration is equal to 

P^EM^E.W^ >i) 

= ^ e tPs 2 Xll}N-s E ^ e t(f3s 2 -N)Xl 1 y E ^-tNXl 1 y(s-l) 

< K[e t^xl l]NE[e -tNXl l]S (s-i) = (1 _ 2tpa 2 )- N '\l + 2tN)-< s - 1 ^ 2 
The right hand side is minimized with 

(l-P)s-l 



t 



2Ps(N + s(s - 1)) 
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and inserting this in the expression 

N + /3s 2 \ ~ N/2 ((s-l)(N + f3s 2 ) \ - s(s - 1)/2 



< 



N + (s-l)sJ \Ps(N + (s - l)s) 

= (N + (s — l) s ) S ( S -l)/2+AT/2 (iv + ^ys(s-l)/2-N/2 ^£__^ ^ ^ 
/ (1-«) S 2\'(-W2 / i\-«(«-l)/2 , w 

< /3 s(s-l)/2 eS /2 e (l-/3)s 2 ( S ( S -l)/2+JV/2)/(iV+/3 S 2 ) < £ S ( S -l)/2 eS /2 e (l-/3)s 2 /2/3 
= e s( S -l)(ln/3)/2+s/2+(l-/3)s 2 /2/3 



If /3 > 1, then we get 

< (1 + 2t/3s 2 )- 7V/2 (l - 2tA0~ s(s ~ 1)/2 
Now the optimal choice for t is 

(l-/3)s-l 
2/3s(iV + s(s - 1)) 

and the result follows 



< 



JV + /3s 2 \ _iV/2 / (s - 1)(N + /3s 2 ) n - s ( s - 1 )/ 2 



^ + (s-l)s/ V^l^ + l 5 -!)«), 

< /3 s(s-l)/2 e s/2 e (l- / 3) S 2 ( S ( S -l)/2+JV/2)/(iV+/3 S 2 ) < ^ S ( S -l)/2 eS /2 e (l-/3) S 2 /2 
= e s( S -l)(ln/3)/2+ S /2+(l-/3)s 2 /2 _ 

Setting i = 1 + 5 or i = + gives the following estimate. 



□ 



Corollary 4.2. Given {Xjj : 1 < j < s, 1 < I < N}, independent identically standard- 
normal distributed random variables, N > s and e > 0, then 

s s N 
VVl 2 ,< 1 C 2S TX 2 ,) < e -s(s-l)(\n(l+8))/2+s/2+6s 2 /2 

3 = 1 1=1 1=1 

< e (l+<5) S /2- S (s-l)(5 2 /2-53/3)/2 

and 

s s AT 

¥(Nj2J2 X h ^ (! + 5 ) s2 E X 5) ^ e »(»- 1 )( ln ( 1 +«))/ 2 +«/ 2 -** 2 /2 

3=1 1=1 1=1 

< e (l+5) S /2- S ( S -l)(«5 2 /2-5V3)/2 _ 

For the next result it is sometimes more convenient to work with random matrices whose 
column vectors are normalized. We will alternate between the random N x M matrix X 
with standard normal entries and the random matrix x which is obtained by 

N 

^ = ^/(E*5) 1/2 
i=i 
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and consequently has column vectors that are independent, uniformly distributed on the 
unit sphere in M. N . 

Theorem 4.3. Let {xi,x 2 , • • • ,x s } be vectors in W N , drawn independently according to a 
uniform distribution on the unit sphere. If V is a fixed subspace of dimension s < N and 
Py is the orthogonal projection onto V , then 

P ( _J_ < Y*\\ PvXi f <1 + 8)>1- 2se^ s ' 2 ~< s - 1 ^ 52 / 2 - 53 ^l 2 . 
l + o ^— ' s 2 
i=i 

Proof. This follows from the fact that mapping a standard normal Gaussian random vector 
to the unit vector in its span leads to a uniform distribution on the sphere. The preceding 
probability estimates are unchanged by scaling the vectors on both sides. Let Zj = Yla=i -^ji 
then the theorem we wish to prove is equivalent to 

P ( _J_ < !L — E Xj,i < (1 + 8)) > 1 - 2s e( 1+<5 ) s / 2 - s ( s - 1 )( 52 / 2 -' 53 / 3 )/ 2 . 
l + o <s Zj ' 

3=1 J 1=1 

To deduce this, we use a union bound again, which implies with the preceding estimates 
that 

1 N s s 

3=1 1=1 

> 1 _ 2se (l+-5)«/2-«(«-l)(<5 2 /2-5 3 /3)/2_ 

Now minj J2i Xf l < Zj < max, J2i X?i establishes the bound. □ 
Next, we prepare the result on the equiangularity of random fusion frames. 

Lemma 4.4. Let M = Ks, M > N > s, < 5 < 1, e = (1 + 5) 3 - 1 and {xj*^ be 
independent random vectors in M. N that are uniformly distributed on the unit sphere. Let 
{Pk\k=\ b e the orthogonal projection onto the span o/{xj}j g j fe , where {</&} partitions the 
index set and each Jk has size s, then for fixed k / I, 

> 1 - R x - R 2 
where the failure probability is bounded by the sum of 

R ± = 2se (l+<5W2-s(s-l)05 2 /2-<5 3 /3)/2 _ 

and 

R 2 = 2(M + K(l + l)-) e -^ 2 /4+^3/ 3 + 2(1 + 4 )JVe _^/4+M*3 /3 



Proof. A random, uniformly distributed s-dimensional subspace is realized by taking the 
span of independent, identically uniformly distributed random unit vectors {xj}| =1 . Corre- 
spondingly, the orthogonal projection is obtained via orthonormalizing these random vectors 
with the square root of the pseudoinverse of S = Ylj=i x j ® x *j- As a consequence, the 
trace is identical to 

s s 

tr[P,P fc ] = J2^[Pl(^) 1/2 Xj x*(S^ 2 } = £ WPiiS^Xjf . 

3=1 3=1 
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With the triangle inequality, we split this expression into three parts that we estimate 
separately, 

Qi - Q2 + Q3 < \\ p i(^) 1/2 Xif <Qi + Q2 + Qs 



with 



and 



Q 3 = EH p '((' st ) 1/2 ^-^)ll 2 - 

The Cauchy-Schwarz inequality gives Q2 < 2(QiQ 3 ) 1/2 so it is enough to control Q\ and 

The quantity Q\ is concentrated near s 2 /N by Theorem 14.31 which gives a lower proba- 
bility bound of 1 - Ri for the set with < < NQi/s 2 < 1 + <5 < 1 + e. The third 
quantity is with probability 1 — R2 small by Theorem 13.41 because if e = (1 + e)) 3 — 1 then 
rr < S < 1 + e and if llxll = 1, then 

\\{S^x - x\\ 2 < (^L= - l) 2 < (vTTe- l) 2 < 6 2 /4. 
yl + e 

Thus, apart from a set of probability given in Theorem El Q 3 < se 2 /A and 2(Q 1 Q 3 ) 1 / 2 < 
(Qisf/h. 

Next, if < Q x < (1 + e)£ then 

1 «2 e 3/2, 

is ,i /oS ' e — „ . ^-^ „ „ , „+si/2 i|2 



ieJfe 
and 

2 3/2 

£ ||PK^) 1/2 ^II 2 < Qi + Q2 + Q 3 < (i + c)^= + (i + ^) 1/2£ t^ + « 2 /4 • 



□ 



By a union bound over all pairings of subspaces, we get the following, more qualitative 
estimate: 

Theorem 4.5. Let K,s,N G N, M = Ks > N > s. For any c > 1 there exists eo > 
suc/i i/iai /or a// < e < eo, if {Vk\k=i ^ s a f am ^V °f s- dimensional subspaces selected 
independently at random according to the unitarily invariant distribution in ~R. N and {Pk\k=i 
the corresponding orthogonal projections, then the set 

[~N N [W Ne 2 

X = {l- ce (l + W_) < ^tr[P fej P,] <l + ce(l + ] J-) + — for all k + 1} 

has probability 

P(X) >l-(Ri + R 2 )K{K - l)/2 
with R\ and R2 as in the preceding lemma. 

Proof. If e is sufficiently small, then the remainder for the series expansion, 1/(1 + e) — 1 + 
e + (VTTe - 1 - e/2)y/N/s is bounded by (c - l)e + ce v / Nj^. A union bound for the 
K(K — l)/2 pairings gives the probability bound. □ 
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Corollary 4.6. Let K, N, s and 5 be as above. If N, s — > oo, s/N — > c > 0, and K is such 
that 

^%^ + ^l„(l + |)<^/4-iV3 

and 

^ln(l + ^)<5 2 /4-<5 3 /3 

then the upper bound R\ + R2 for the probability of X in the preceding theorem decays 
exponentially in N . 

Proof. As s/N remains bounded away from zero as N — > 00, Ri decays exponentially. 
Examining the exponents of the terms in R2 shows that if the conditions on K, s and N 
are satisfied then it decays exponentially as N — > 00 as well. □ 
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